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Abstract 

o 

' We propose a novel basis of vector functions, the mixed vector spherical harmonics that are closely related to the 
functions F\ m of Sheppard and Torok and help us reduce the concentration problem of tangential vector fields within 
a spherical cap to an equivalent scalar problem. Exploiting an analogy with previous results published by Griinbaum 
, and his colleagues, we construct a differential operator that commutes with the concentration operator of this scalar 
problem and propose a stable and convenient method to obtain its eigenfunctions. Having obtained the scalar eigen- 
£\j , functions, the calculation of tangential vector Slepian functions is straightforward. 
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1. Introduction 

In our mathematical models, we often assume bandlimitedness for physical or computational reasons, yet also wish 
that our solutions be localized, with respect to their energy, inside a finite spatial region. Since these are mutually 
exclusive conditions |l|], we need to resort to bandlimited functions with an optimal spatial localization. The goal of 
the so-called spatial concentration problem Q P- 75] is to find such functions, and since its first thorough investigation 
\Q ' by Slepian, Landau and Pollak in one and multiple Cartesian dimensions Orel], it has been revisited many times, 
CNI , including solutions for spherical and planar regions of arbitrary shape ll7l-[loll. 

Each individual concentration problem gives rise to an orthogonal set of well-localized functions, which now we 
refer to under the common name Slepian functions. They have enjoyed increasing popularity in applications involving 
signal processing, function representation and approximation or the solution of inverse problems. In particular, the 
scalar spherical Slepian functions have been utilized, for instance, in geodesy and geophysics JI - 12 ] , cosmology 1 13 
[Till , computer science lfl5ll and mathematics lfl6ll . 



While they have been widely applied in the last two decades, it was not until recently that the theory of vector 
Slepian functions began to mature. The first successful construction of bandlimited vector fields, localized to a spher- 
ical cap, was reported in the context of biomedical science lfl7ll . followed by an application in physical optics lfl8ll . 
Recently, a more general treatment of the vector spherical concentration problem has been published [ 19], however, 
the question on the existence of a commuting differential operator for the spherical cap has been left unresolved. This 
question is important for the following reasons. 

Slepian functions of a particular problem are eigenfunctions of the so-called concentration operator associated 
with the spatial region of interest, an integral operator exhibiting a peculiar step-like eigenvalue spectrum. This 
property makes the direct calculation of its eigenfunctions numerically difficult B20I1 . 

A particularly important result of Slepian and his colleagues was the introduction of a Sturm-Liouville differential 
operator that commutes with the corresponding concentration operator, hence they share a common set of eigenfunc- 
tions. Since the differential operator has a simple spectrum with more evenly distributed eigenvalues, it allows a 
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more stable and accurate numerical computation of the eigenfunctions |20]. Two decades after the seminal papers 
by Slepian and his colleagues, Griinbaum, Longhi and Perlstadt found such a commuting differential operator for the 
scalar concentration problem within a spherical cap as well 17:]. However, we are not aware of a similar proposal for 
the vectorial problem. 

In this paper, we construct a differential operator commuting with the concentration operator of the vector case. 
This constitutes our main result. We restrict ourselves to tangential vector fields, since the concentration problem of the 
radial component is equivalent to the scalar concentration problem on the sphere II 1911 . which has already been studied 
extensively [9]. The key functions in our investigations are the novel mixed vector spherical harmonics Qf ln (6,cf>) 
which enable us to reduce the vectorial problem to a scalar one involving the special functions F[ m of Sheppard and 
Torok 12111 . After that, the problem can be solved in an analogous way to its scalar counterpart |9J. 

We note that we only consider spatially localized, bandlimited fields here. The symmetric problem of spectrally 
concentrated, spacelimited functions can be derived by exploiting the analogy to previously published results fl9ll . 



2. Preliminaries: associated special functions 

In preparation for the concentration problem, we give a detailed survey on the essential properties of important 
special functions used in our investigations. We start by revisiting the key properties of the normalized associated 
Legendre functions, because they provide a foundation for the theory of the special functions F/ m . After that, we 
prove several fundamental relations for F/ m , such as orthonormality and recurrence relations, and Christoffel-Darboux 
formula. 

Next we diagonalize the vector Laplacian on the spherical surface and introduce the mixed vector spherical har- 
monics. We establish a relation between them and the functions F/ m and discuss their orthogonality properties. Finally 
we show how to expand an arbitrary tangential vector field in terms of the mixed vector spherical harmonics and define 
bandlimitedness in this context, so that we can use these new vector fields as basis functions in the treatment of the 
concentration problem in Section [3] 

2.1. Normalized associated Legendre functions 
2.1.1. Definition and orthonormality 



The normalized associated Legendre functions Ui,„ of integer degree / and order m are defined as [22J, p. 757] 



U lm (x) := c lm P'l'(x), I > 0, -Z < m < I, (1) 

where 



(— 1 \ m c\ ,+m 

( -iL ( i-x 2 )^ 2 -, 

2 l l\ dx l+m 
are the (unnormalized) associated Legendre functions [22, p. 743], and 



Pf( x ) := -^-(1 - x 2 ) m/2 ^ir(* 2 - 1)' (2) 



2/+l(/-m)! 



\ 2 (l + m)\ 



is the normalization factor. 

The normalized associated Legendre functions satisfy the orthonormality relation 



X 



1 

Ui m (x)Uv m (x) Ax = S u >, (4) 



where 5w is the Kronecker delta. 
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2.7.2. Recurrence relations 

We can obtain recurrence relations for U[ m in a straightforward way by transforming the corresponding relations 
for P™ |22l p. 7441: 



xUi m (x) = & m Ui-i >m (x) + &+\, m Ui +Um (x), 

= -lxUi m (x) + (21 + Y)t;i m Ui-i m (x), 



(1 -x ) — 

ax 

2 dlli m (x) 

(I -x ) — — — = (/ + l)xUi,„(x) - (21 + l)£i +hm Ui +hm (x), 



where 



Ax 



&m '■ = 



I + m ci„ 



(I + m)(l - m) 



2/ + lQ-i,,„ \(2/+l)(2/-l) 



In addition, two more recurrence relations can be formulated II23L 12411 : 

^-x 2 U '!" ^ = al m U Lm+1 (x) + aJ m U l x (x), 



Ax 



where 



and 



where 



y/(l + m)(l + m + 1) 



mUi, n (x) 
Vl -x 1 



= b tm U l-l,m + l( x ) + b lm U l-\M-l^ 



lm 



21+ 1 ^(l + m)(l + m- 1) 



(5) 
(6) 

(7) 
(8) 

(9a) 

(9b) 
(10a) 

(10b) 



2/- 1 2 

Note that everywhere in this paper, where the ± or + signs occur, either the upper or the lower one has to be used 
consistently in the whole expression. 

2.1.3. Differential equation and symmetries 

The functions [//„, are solutions to the Sturm-Liouville differential equation 11221 p. 744] 



d 

dx 



(i _ x2) d«(x) 
dx 



1 — x : 



■u(x) — —1(1 + l)u(x), 



(11) 



known as the associated Legendre equation. Since it contains m 2 only, t// m and Ui- m must be proportional. In fact, 
they are related by the symmetry relation I122[ p. 743] 



Ui,- m (x) = (-\) m U lm (x). 



In addition, we can formulate the parity relation [22, p. 746] 

U, m (-x) = (-l)' +m U Im (x). 



(12) 



(13) 



2.1.4. Special values 

There exist closed-form expressions for special arguments or parameter values of Ui,„. At the interval endpoints, 
Ui m evaluates to 1221 p. 746] 

f(±l)'c y , ifm = 0, 
1 otherwise. 



U,,„(±l) 



(14) 



Another expression for the case I — m is [22, p. 745] 

U mm (x) = (-l) m c mm (2m - 1)!! (1 - x 2 )" ,/2 , m > 0, 
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(15) 



where (2m - 1)!! = (2m - l)(2m - 3) • • ■ (1) is the double factorial. Equation (T5[ can be used as an initial value 

i I 

to obtain Ui m (x) numerically in a fast and stable way 11251 p. 364]. Considering m > 0, we set C/ m -i, m (x) = and 
calculate U mm (x). We then use recurrence relation (0 in the upward direction until we reach Ui m (x). Function values 
for m < can be obtained using symmetry relation (TTZV 

2.1.5. Addition theorems 



Based on the work of Winch and Roberts 126J, we can formulate two addition theorems which will be useful later: 



dU,,„(x) 



dx 



^ [mUi m (x)] 



l(l+\)(2l+ 1) 



2 ;<7+l)(2/+l) 



(16a) 
(16b) 



m=-l 



2.2. Special functions Fi m 

2.2.1. Definition and orthonormality 

In this section, we give a detailed description of the functions 



Fim(x) := 



1 



VZ(TTT) 



Vl -x 2 



dUl m (x) 

dx 



vr 



:Ui m (x) 



I > 1 , —m < I < m, 



(17) 



previously defined by Sheppard and Torok [21]. We can express the conditions for the integer indices I and m alter- 
natively as -oo < m < oo and / > t m , where the minimal degree ( m for a fixed m is 



£ m :- max(l, |m|). 
The functions F\ m are orthonormal over [—1,1] for fixed m, i.e. 



x: 



Fbn(x)F Vm (x) dx = 5 W , 



(18) 



(19) 



as shown in |Appendix A. 1 1 A small subset of them (I < 3) is described in Fig.Q] 

We can obtain useful equivalent formulations of F; m by inserting recurrence relations (0 or ([7} of U\ m into (TTTl : 



Flm(x) 
Flm(x) 



-(lx + m)Ui m (x) + (21 + l)& m Ui-i, m (x) 
V/(/+ l) Vl -x 2 

[(I + \)x - m]U tm (x) - (21 + l)^ /+ i, „,£//+! ; „,(x) 

V/(T+ l) Vl -x 2 



(20) 
(21) 



Sometimes it is inconvenient that expressions (fTTT i. d20l i. and OTb are all singular at x = ±1 because of the factor 
(1 - x 2 y l l 2 . However, a singularity-free form can also be obtained by exploiting recurrence relations (|9} and dT(Tt : 



Fl,„(x) = 



-1 



V7(7TT) 



(22) 



It is straightforward to show using (TTZt and (l22t that for the special case m = 0, the equivalence Fifl(x) = L^i(x) 
holds. 
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2.2.2. Recurrence relations and Christoffel-Darboux formula 

We have found the following recurrence relations for Fi„, (cf. the corresponding relations <(5)-(|7} of £//,„): 



1(1+ 1) 
(1-x 2 ) 



Flm(x) - £l m Fl-\jn( X ) + (l+\,mFl+\, m (x), 

dF tm (x) 



dx 



, dFi m (x) 
(1-x 2 ) — l f^ = (l+l) 
dx 



X ~~ p. ) F l™&) + ( 2/ + l)£lmFl-l, m (x), 

Flm(x) - (21 + l)(M, m Fi +lm (x), 



(/+D 2 



where 



_ v(/+i)(z-i) ^ _ i 



(/+ 1)(/- l)(/ + m)(/-m) 



(2/+ 1)(2Z-1) 



(23) 
(24) 
(25) 

(26) 



A proof for (1231 and (1241 is provided in Appendix A. 2 and Appendix A. 3 respectively. Relation (l25l is straightfor- 
ward to prove by combining (1231 and (l24l . 

Relation d23l can also be used to derive a Christoffel-Darboux formula fl27l p. 42] specialized for F; m : 



(X - x') 2^ F, m (x)Fi m (x) = (L+l,m[FL+l,m(x)F Lm (x') - F Lm (x)F L+hm (x')]. 



(27) 



See Appendix A. 4 for the details. 

2.2.3. Differential equation and symmetry 

The functions F/ m satisfy the Sturm-Liouville differential equation 



d 

dx 



2 du(x) 
(\-x )— — 
dx 



m 2 - 2mx + 1 
1-x 2 



-m(x) = — /(/ + I)m(x), 



(28) 



as proven in Appendix A. 5 Comparing its symmetry properties to those of the associated Legendre equation (fTTb . 
we find an important difference: while (fTTT > is invariant under a change in the sign of m, in (l28l > the signs of both 
m and x have to be changed for transformation invariance. In fact, by inserting symmetry relation (1T21 of £//,„ into 
definition ( fTTT i of F/ m and exploiting the parity relation ([T3l . we get the following symmetry relation for F; m : 



F,,- m (-x) = (-l)' +1 F, m (x). 
This symmetry relation is also apparent in Fig. [1] 



(29) 



2.2.4. Special values 

Combining i// m (±l) of (TT4l and the singularity-free form (l22l of F/ m provides a simple way to calculate the 
function values at the endpoints x = ±1: 



F fa (l) = 
Fm(-1) = 







if m = 1 , 
otherwise, 



(-D'-'q. 




o 



if m — - 1 , 
otherwise. 



(30a) 
(30b) 



In addition, to get a relation similar to the closed-form expression (IT5b for t/ mm , we can combine ( fT3b with 
definition ( fTTb of F/ m . In this way we obtain 



F( m ,m(x) 



(-i r +i (i + x)(D m (x) 

(1 -x)<D w (x) 
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if m > 0, 
if m — 0, 
if m < 0, 



(31) 



where 



<S> m (x) :-- 



171+1 



c mm (2m+\)\\{\-x i ) 



^(m-l)/2 



(32) 



Like in the case of the associated Legendre functions, recurrence relation (1231 provides a stable and efficient 
method to evaluate Fi m (x) numerically. By setting Fi m -\ m = and starting with the closed-form expression (|3H for 
Fe mt m(x), recurrence relation (l2Jt can be used repeatedly in the upward direction until one obtains Fi m (x). 



2.2.5. Addition theorem 

We can use the addition theorems ( TToT l of Ui m (x) to formulate a similar relation for Fi„,(x) as well: 



[F lm (x)f = 

m-—l 



21+ 1 



(33) 



The proof can be found in Appendix A. 6 



2.3. Scalar and vector Laplacian on the unit sphere 

Let // = u(9, <p) be an arbitrary scalar field and v = Vg(6, <p)6 + v$(6, 
unit sphere 

Q := {(9, (p) : < 9 < n, < (p < 2n) . 

Here, 6 and are unit vectors in the 6- and ^-directions, respectively, as shown in Fig. [2] 
The surface scalar and vector Laplacian on Q are defined as follows [28]: 

A2' 



a (tangential) vector field defined on the 

(34) 



1 d 



sin 9 89 



8u , 
sin 9— + 

89 



1 



sin 2 



cos <9v0 



sin" 



sin' 



dcp 



+ 



+ 2 



cos 9 dve 



sin 



2 9 sin 2 50 



(35) 
(36) 



Note that the each component of V 2 ^ v contains both components of v. We can diagonalize V 2 ^ v by introducing the 
tangential basis vectors 



1 

V! 



+ i4>), 



which are orthogonal with respect to the complex dot product 

t : ■ t x = o. 



(37) 



(38) 



Here i is the imaginary unit and the asterisk denotes the complex conjugate. 
In this new basis, we can write (l36T l as 

V 2 1 v = (A+v + )T + +(A n v_)t-, 



(39) 



where v = v+(6, <p)f+ + \'-{9, <p)f~ and 



A* •- V 2 - 



sin 2 



1 ± 2icos#— 



(40) 



Next we introduce two fixed-order operators, which play a central role in constructing commuting differential op- 
erators for the spherical cap, as we will see in Section |3~3| For separable functions of the form u(9, (f>) - w(9) exp(im^) 
the differentiation with respect to <p can be performed explicitly, and the scalar operators V 2 ^, A+ and A n become 
identical to the following fixed-order operators V 2 , , A nm and A n , respectively: 



1 d 
sin 9 d9 



A •- V 2 



sin6»— I - 
d6j sin" 

1 - 2(+m) cos 9 
sin 2 9 



(41) 
(42) 



Since V 2 , = _ , we can write ( l42l in a more compact form: 



A •- V 2 - 



1 - 2m cos 



(43) 



Upon substituting x = cos 6 in V?, and A n , we regain the differential operators on the left-hand side of differ- 



ential equations ( fTTb and ( l28l i. respectively. Therefore £//,„ and Fi m satisfy the eigenvalue equations 

V^ m t/, m (cos0) = -/(/+ l)I/ /m (cos0), 

A n , m .F/ m (COs6>) = -/(/ + l)F/ m (C0s6»). 

2.4. Mixed vector spherical harmonics 
2.4.1. Definition and orthonormality 

We define the mixed vector spherical harmonics as 

Qf m (6, 0) := ^^[Y /m (0, <p) ± i Z /m (0, 0)], 
V2 



(44) 
(45) 



(46) 



where Y/ m and Z/ m are the conventional (fully normalized) tangential vector spherical harmonics. They are defined 

asESB 



Yi m (6, </>):= 
Z lm (6,4>):= 



VZ(TTT) 
i 



i dY lm (9,4>) h dY lm (e,<p)~ 

-" — <v 



sin 6 d<f> 



06 



dYUe,4>) h 1 dY lm (0,4>). 

u + 



where 



y/(z + 1) 

F /m (0,0):=t/ /m (cos0) 



56* " ' sine 
exp(im^) 



V2^ 



(47a) 
(47b) 

(48) 



are the scalar spherical harmonics. 

Equations ( TTTl i. ( l37l >. ( l47i i and (l48b can be used to write (PES! in a separable form: 



Qfa(M) = /W«w$> 



exp(im</ 
V2^ 



(49) 



It directly follows from this formulation and the orthogonality of f ± that, unlike Y/ m and Z/ m , the functions Q^ n exhibit 
local (vector) orthogonality, regardless of their degree and order, i.e. 



Q*(M)-<&,(M) = 0. 

This equation together with (l49l can be used to prove the orthonormality relations 



r 



Qf m *(0,0)-Q* m ,(0,0) dn = 6v6 K 



(50) 

(51a) 
(51b) 



where J^... dO := . . . sin6»d6>d0 
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2.4.2. Special values 

The values of Qf m at the ^-coordinate singularities deserve extra attention, since f± are not well defined there. 

We can circumvent this problem by expressing Qf using Cartesian basis vectors (see Fig. Using Eqs. (f30T > for 
F, m (±l), we have 

I t^=c/o(x ± if) ifm = ±l, 
Qte(0 = O,0) = L . (52a) 

10 otherwise, 

f^pQo(x + iy) if m = +1, 
Qf m (6 = *,<[>) = ' ■ ( 52b ) 

10 otherwise. 



2.4.3. Spherical harmonic expansion and bandlimited functions 

Like Y/,„ and Z/ m ll30ll . the mixed vector spherical harmonics also form a complete basis of the Hilbert space 
of square-integrable tangential vector fields defined over Q. Hence we can expand an arbitrary tangential vector field 
v in terms of Qf m as 

00 / 

W <(>) = 2 2 [vLQtje, <f>) + vJMe, </>)] , (53) 

1=1 m=-l 

where the expansion coefficients vf m can be calculated as 

vf m -= f Q£T(0,0-v(0,0)dn. (54) 
Jn 

If vf m = for L < I < oo and some L > 0, we call v bandlimited. The limit L is the maximal degree of functions 
Q^„ that contribute to the expansion (f53t . Therefore the subspace S l of bandlimited vectorfields is finite dimensional, 
its dimension is equal to the number of terms in (l53l : 

L 

dimS L = 2 YP- 1 + !) = 2 [( L + !) 2 ~ !] = 2L ( L + 2 )- (55) 

l=i 



3. The concentration problem of tangential vector fields within an axisymmetrical spherical cap 

Having defined all important special functions, we now turn our attention to the main topic of the paper. After 
a brief discussion of the concentration problem in terms of the mixed vector spherical harmonics, we introduce a 
reduced scalar problem which can be solved analogously to the theory of scalar spherical Slepian functions 17l|9|]. 

In Section [3721 we analyze the eigenvalue spectrum of the concentration operator and give an illustration on the 
scalar eigenfunctions. 

Finally, we propose a fast and numerically stable way to calculate the eigenfunctions by using a differential 
operator that commutes with the scalar concentration operator obtained previously. 

3.1. Formulation of the vector problem and its reduction to a scalar one 

We consider the variational problem of finding a bandlimited, tangential vector field G{0, <p) e S l that maximizes 
the fractional energy contained within an axisymmetric spherical cap C: 

f,|G(0,0)| 2 dQ JT G*(0,0)-G(0,0)dfl 

max — = max — . (56) 

g J n \G(0,<p)\ 2 dn g J n G*(0,0)-G(0,0)dQ 

Without loss of generality, we can center our spherical cap at — 0, as seen in Fig. (|3): 

C = {(0,4>) : < 9< 0,0 < (f> < 27!-), (57) 

where > is assumed. 
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To solve the vectorial problem ((56), we adapt the method of the scalar case Js}] and turn ((56) into a Rayleigh-Ritz 
matrix variational problem 113 ll p. 176]. We can achieve this by expanding G in terms of Q^ n : 



L I 



G(6, <j>) = J] Z [stnQtn^ 0) + giStiJA = Z Z [gtmQUO, 4>) + SfaQFrn^ 0)] • (58) 



1=1 m=-l 



In the second step, we have interchanged the order of double summation to facilitate the transition to the matrix 
formulation of (156) later. A visual comparison of the two summation schemes is given in Fig.|4] 

Next we insert (158) into ((56), interchange the order of summation and integration and use the orthonormality 
relations (150) and (IBD . Hence the numerator can be written as 



6 I'm' 



XLt Li Li Li V s~\ 

G*(e,0)-G(0,0)do= 2 Z^* Z Z Q£0.0- 0<Jfl 
m=-Ll=e m m'=-Ll'=e m , L^ c 

L L L L r n 

+ Z Z^» Z Z Q^e^-Qj^ie^dQ g j, m „ 

m=-Ll=t,„ m'=-L !'=(,„, *- Jc 

while the denominator becomes 

XL L 
G\6, cf>) • G(6, </>) dQ = J] J] + . 

The integrals on the right-hand side of (|59) can be expressed as 

J Qte (0, 0) ■ QL'( - 0) d ^ = ^ J[ exp[i(m' - m)<£\ d<p jT F,, ±m (cos 6»)F P , ±m (cos 0) sin 0d0. 

Since (27r) _I C 2 * exp[i(m' - m)0] d</> = (5 mm <, Eq. (loD further simplifies to 



(59) 



where 



A", 



Jo 



F fa (cos 6)Fi, m (cos 6) sin d0 = I Fi m {x)F Vm {x) dx. 

J cos© 



(60) 



(61) 



(62) 



(63) 



The integrand of K m ji' is a polynomial of degree I + I', hence it can exactly be integrated numerically, for instance, by 
a Gauss-Legendre formula of [{I + I' + l)/2] nodes. 
Taking d62) into account, (156) can be rewritten as 



max 

total 



l ( L 



l ( L 



Z Z Z ^.«'^ + Z Z g to Z K -m,U'gJ'm / Z Z (^toT^toi + S/mSfa) 



m=-L\l=( m I'-. 



m=-L\l=t m I': 



m=-L !=(,„ 



(64) 



To express (164) in matrix formalism, we construct a column vector g of 2L(L + 2) elements as 

9 := \gL,-L>gL-l,-L+l>gL,-L+l> ■ ■ ■ >gL,L'gL,-L>gL-l,-L,+l>gL,-L,+l> • • ■ >8l,L\ ■ 

where the expansion coefficients gf are enumerated according to the scheme of Fig. 0?. In addition, we introduce the 
2L(L + 2) x 2L{L + 2) block-diagonal matrix 

K+ 
K 



K 



(65) 



(66) 



where the L(L + 2) x L(L + 2) blocks K ± are themselves block-diagonal: 

K + := diag [K_ L ; K_ i+1 ; . . . ; K L ; ] , K~ := diag [K L ; K w ; . . . ; K_ L ] . 
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(67) 



The elementary blocks 

K m i i ■ • ■ K m 



K,„ - 



K m ,LC w ' ' ' Km,LL 



(68) 



correspond to different orders to and have an order-dependent size of (L - i m + 1) x (L - ( m + 1). Note that order of 
the blocks K m in K ± is reversed owing to (l62l l. 

We can thus use these constructions to transform problem (l64l into the Rayleigh-Ritz matrix variational problem 

max (69) 

9 g 1 ' g 

where the dagger sign denotes the conjugate transpose. Equivalently, we have to find the eigenvector g of the eigen- 
value problem Oil p. 176] 

Kg = 779 (70) 

with the maximal eigenvalue 77. However, rather than solving the large 2L(L + 2) x 2L(L + 2) eigenvalue problem (l70l . 
the block-diagonal structure of K allows us to solve a series of smaller (L — € m + 1) X (L — l m + 1) problems instead, 

K m g m = VmQm, ~L < m < L, (71) 

one for each order m. 

From d63j follows that K m ji> = K m jq, implying that the matrices K„, are symmetric. Hence their eigenvalues are 
always real. For a given to, we rank-order the (L - l m + 1) distinct eigenvalues r\ mn as 1 > r} m> \ > t] m ^ > ■ ■ ■ > 
T)m,L-t„+i > 0. The associated eigenvectors g,„„ can be chosen to be real and orthonormal: 

9ln9mn> = 5 nn>> OLKiQmn' = 1mAn>, 1 ^ ». n' < (L - € m + 1). (72) 

Here we have distinguished between the different eigenvalues and the corresponding eigenvectors by the use of the 
additional index n (or «')■ However, we drop this additional index for brevity when we refer to any of the (L — l m + 1) 
eigenvalues or eigenvectors. 

We also denote the elements of an eigenvector g,„ simply by gi,„. This brings up the question: how are the 
coefficients gi m connected to the original coefficients gt and gj of expansion ( T58b ? According to (|66"I > and d67l i. each 
block K m occurs twice in K, hence each eigenvector g m gives rise to two vectorial eigenfunctions: 

L L 

G + je, <P) = Yj SIQIJO, <p) = J] gbn Qt,(e, <f>), (73a) 

L L 
GZ m {6, 0)=2 Sl-n&l-M *) = E SlmQJ-mi ' 4>)- (73b) 

Its worth emphasizing that every eigenfunction contains either Q^, or Qj m of a single order to only, which is a 
consequence of the block-diagonal nature of the concentration matrix K. Upon substituting expression (l49l of Q^ J; 
into Eqs. (1731 , we obtain 

G^(M)=CWcoB0^^t ± . (74) 

where 

L 

G,„(x) := ^ gimFimix) (75) 

are real functions. 
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In this way, we managed to reduce the vectorial concentration problem within a spherical cap to equivalent one- 
dimensional, scalar concentration problems of various orders m. The key idea in this simplification was the choice d49l 
for our basis functions. The scalar concentration problem for a fixed order m can be formulated as 



max ■ 



£ osQ [GJx)fdx ^ 
f^[G m (x)] 2 dx ' 



(76) 



where G m is a bandlimited scalar function belonging to the subspace spanned by Fi,„. The corresponding Rayleigh- 
Ritz matrix variational problem is 

«T is « 

y*;.' infill (77) 



max ■ 



9'" 9m9m 

Instead of the eigenvalue equation ( 1711 ) specifying eigenvectors g,„, we can directly formulate an eigenvalue equa- 
tion in terms of the functions G m , too. Therefore we first express Eq. (fTTl component-wise as 



/ \ K m ji'gi> m - rj m gi m , l m <l < L. 



Now we multiply both sides by Fi m (x) and sum over I: 

L L L 
/ \ / l K m Jl'gl> m Fi m (x) = T] m glmFl„,(x). 



(78) 



(79) 



i=t,„ i'=t, 

The left-hand side can be rewritten as 

L L 



l=tm 



L L L L v ~\ 

^ 2 K m jvgv m Fi m (x) = ^ ^ I FiJx')F hn (x') Ax' 
l=t m l'=t m L-'™ 50 



gl'mFlmix) 



J J t 



, Fi m (x)F tm (x') 



2^ gl'mFl'm(x')[dx'. 



l'={,„ 



This way, we obtain a Fredholm integral equation of the second kind for G,„, 

% n (x,x')G m (x')dx' = rj m G m (x), -1 < x < 1, 

) 

where the kernel function < K m is defined as 



Jcos© 



% n (x,x') := ^ F tm (x)Fi m (x'). 



(80) 



(81) 



It follows from the orthogonality relations ((72} of the eigenvectors that the scalar eigenfunctions G m are doubly 
orthogonal: 



I G mn ^X)G mn \x)dx — 6 nn ', 

J-l 

f 1 

Jcos & 

The vectorial eigenfunctions G* inherit this property as well: 

f G±„*(0, 0) • G% n ,{9, cf>) dQ. = <W<W, f G±„*(0, 0) • G*,„,(0, 0) dO = 0, 
Jn Jn 

£ G^(0, 0) • G% n ,{6, cf>) da = /7™<W<W, G^(0, 0) ■ G* „,(0, 4>) dQ. = 0. 



(82a) 
(82b) 

(83a) 
(83b) 
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3.2. The eigenvalue spectrum and its peculiarity 

The eigenvalue spectrum of Slepian-type concentration problems JUSll^l exhibits a characteristic step-like shape, 
and the present case is no exception. Figure [5] shows rank-ordered spectra including rj mn for all orders m. They 
correspond to = 30°, 60°, 90° and the maximal degree was chosen L — 18. 

The majority of the eigenvalues for each case is either close to one or zero, corresponding to well-concentrated 
and poorly concentrated eigenfunctions, respectively. As an illustration, in Fig.|6]we have plotted a small number of 
scalar eigenfunctions G mn , corresponding to different parts of the eigenvalue spectrum. 

Strictly speaking, the solution of the concentration problem (l56l > is the pair of vectorial eigenfunctions which 
corresponds to the maximally concentrated G,„. However, having solved the equivalent eigenvalue problem (l70l . we 
have gained a whole set of well-concentrated, orthogonal pairs of eigenfunctions G* . How many pairs belong to this 
set? To answer this question, we first define the partial Shannon number |9[] 

N m :=TtK m = V rjmn = I e K m (x,x)dx, (84) 



which gives the approximate number of reasonably well-concentrated (r] > 0.5) scalar eigenfunctions for a given 
maximal degree L and order m. Summing over all possible values of m, we obtain the (total) Shannon number 

l l L-e„,+\ _i l ^ 

N:=YjU m =Yj Yj n*m= \ Y J 'K m (x,x)dx = L(L + 2)^, (85) 

m=-L m=-L n=l ^ cm@ m=-L 

where Ac = 2n{\ - cos©) is the area of the spherical cap C. In the last equality we substituted definition ( |8TT >. 
interchanged the order of double summation and used addition theorem ( 133) . 

Hence there are N pairs of orthogonal vectorial eigenfunctions which are suitable for approximating bandlimited, 
tangential vector fields localized to C. Equivalently, the use of this basis reduces the number of degrees of freedom 
from dimS L = 2L(L + 2) to 2N. 

3.3. Toward an efficient numerical solution: the commuting differential operator and its eigenvalue problem 

In Section [3~T1 we obtained the expansion coefficients gi„, by solving eigenvalue equation (TTTb directly. However, 
while it is theoretically possible to calculate gi„, this way, the accumulation of the eigenvalues rj at one and zero, as 
seen in Fig. [5] makes the numerical solution of ( TTTT > ill-conditioned [20]. In order to circumvent this problem, we set 
out to construct another matrix with a simple spectrum to supply the expansion coefficients g\ m . 

Therefore, we first return to the Fredholm eigenvalue equation (l80l . We wish to find a Sturm-Liouville differential 
operator J m that commutes with the concentration (integral) operator on the left-hand side of (l80l : 

f % n (x,x')J' n u(x')dx' =J m f <K m {x,x')u{x!)&x' = f J m % n (x,x')u(x')dx' (86) 

JcosQ «7cos© «7cos@ 



for any square-integrable bandlimited function u, so that the two operators share a common set of eigenfunctions fl22 



pp. 314]. It is known from the Sturm-Liouville theory that J" m has a simple spectrum of distinct eigenvalues with an 



accumulation point in infinity B32l p. 724]. If such a differential operator J~ m can be found, its matrix representation 
can be used to obtain the expansion coefficients g/ m (hence the eigenfunctions) in a numerically stable way. 

The same approach was taken by Griinbaum et al. for the concentration problem of scalar functions within C J3l. 
They proposed the differential operator 

, d 

Q m : = (cos © - cos 0) Vl m + sin 6 — - L(L + 2) cos 6, (87) 

da 

where m is the fixed-order surface scalar Laplacian (HTt . This operator commutes with the concentration operator 
of the scalar case which contains the kernel function T) m {x, x') - Y!{ = \ m \ Ui m (x)Ui m (x') J9J]. 
Based on d87b . we make the following ansatz on $ m ; 

d 

J_ := (cos - cos 6) A n m + sin L(L + 2) cos 6, (88) 

d6> 
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where An m is the fixed-order operator (l43l related to the surface vector Laplacian over SI. Changing the variable to 
x = cos 8 yields 



J m = (cos - x)An,m -(}- x )—-L(L + 2)x, 

ax 



which is equivalent to 



d 

dx 



(cos©- jc)(1 -X-) 



dx 



L(L + 2)x - (cos — x)- 



m — 2mx + 1 



1 -x 2 



(89) 



(90) 



To prove that satisfies the commutation relation (l86l l. we suggest following the concept of Griinbaum et al. J3l. 
First, one proves the identity 



I M] 

J cos© 



(x) [J m u 2 (x)] dx 



= I [Jm«1 
J cos 



(X)] Mafjc) cbc, 



(91) 



which holds for any two functions u\ and u 2 that are non-singular at the interval endpoints (see Appendix A. 7 for 
details). Therefore, the left-hand side of the commutation relation (l86T l can be rewritten as 



f % n (x,x')J' m u(x')dx' = f [X7U^')]«(*')d*'. 

Jcos© J COS® 



(92) 



Finally, one verifies that 

J,Mx, x') = X^m(^. *')■ (93) 

The proof of d93l . like the proof of (|9TV closely resembles its counterpart from the scalar concentration problem [0. 
The key steps are the same, with the main difference that the associated Legendre functions are replaced by F/ m 
together with the corresponding identities. The details can be found in |Appendix A.8| 

Since J',,, commutes with the integral operator of (l80l , the functions G m are eigenfunctions of JT m , too: 



J„,G m (x) =XmG m (x). 



(94) 



Figure Q shows the ^-eigenvalue spectrum of all orders m for = 30°, 60°, 90° and L — 18 (cf. Fig. [5). Similarly 
to the scalar concentration problems (jj, the rank-ordering for Xmn is the opposite of the rank-ordering for rj mn . 
Importantly, the ^-spectrum does not exhibit an accumulation of eigenvalues. 

To obtain a matrix equation similar to the component-wise eigenvalue equation (1781 of K m , we substitute expan- 
sion (fTBl of G\ m in terms of F\ m into eigenvalue equation (1941 . but this time, writing I' instead of /. After that we 
multiply by F\ m (x), integrate over -1 < x < 1, and invoke orthonormality relation (fT9t of F\ m . In this way, we arrive 
at the equation 



where 



v=t m 



J, 



,,«':= J F lm (x)J m F Vm (x)dx. 



Similarly to K m , we can arrange J m ji> into a matrix J m : 



Mm 



J n 



(95) 



(96) 



(97) 



However, the only non-zero matrix elements, as proven in |Appendix A. 9 

J m,ll 



are 



-/(/ + l)cos@ + m 



1 



L(L + 2) + 1 



1(1+1) 

Jm,l,l+\ = Jm,l+U = [1(1 + 2) - L(L + 2)]£/ + i jm , 



(98a) 
(98b) 
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hence J m is real, symmetric and tridiagonal. The eigenvalue equations d95l ) can thus be written as 

*}m9m=Xm9m, -L <m < L. (99) 

We have already seen in Section l2.2.1l that F/o - £/;,i, hence in the special case of m = 0, matrix Jo is identical to the 
matrix of the Griinbaum operator Q\ J9J]. 

In summary, to calculate the scalar eigenfunctions G m for each order m, we first construct the tridiagonal matrices 
J,„ using formulae ((98) and then solve the corresponding eigenvalue problem (|99l numerically. The resulting eigen- 
vectors g m contain the expansion coefficients gi m , i m < I < L, which, substituted into expansion d75l l give the eigen- 
functions G m . The corresponding energy concentration ratio rj m can be calculated using either r\ m — l zo ^QiG m {x)^ Ax 
on] m = g^,K m g m . 

Finally, we demonstrate the numerical stability of the proposed method. We calculated the eigenvectors Q m „ for 
m — 1, L = 18 and © = 30°, 60°, 90° in multiple ways. First, as a reference, we used arbitrary precision arithmetic to 
obtain the eigenvectors of Ki with the relative error of each coefficient gi \ being less than 10~ 23 . Let denote these 
vectors. Then we computed both Kj and Ji in double precision and fed them into the divide-and-conquer routines of 
LAPACK 1 33] to produce the eigenvectors again. Let and g] stand for these results, respectively. In addition, 

we furthermore assume Hg^ll = llgfjl = JI = 1 where ||v|| := Vv T v. 
Figures HJa-c) plot the eigenvalue gaps [33j p. 104] 

gap(77 1>n ) := min|77 1>n - 771 J, (100a) 

j*" 

g a PCti,«) : = min \XU -XiJ, (100b) 
for all three values of 0, respectively, where 1 < n < L. Figures, ©d-e) show the errors 

*g?„ := min(||g^ - g£||, ||g^ - (-g^)||) , (101a) 
6qI := min(||g J Ul - g^||, ||g J ljf; - (-gf„)||) (101b) 

of the eigenvectors, where we have taken their sign ambiguity into account. 

In Figs.Oa-c), we clearly see the accumulation of eigenvalues tjij, of Ki for both small and large values of n. The 
decrease in the eigenvalue gap by many orders of magnitude is accompanied by a rapid increase in the error 6gf lf33l 
p. 104], as seen in Figs.[8jd-f). Therefore, with a naive treatment of Ki, we failed to calculate the well-concentrated 
eigenfunctions accurately; precisely those that play an important role in the approximation of functions localized to 
C. 

On the contrary, Figs. |8]a-c) demonstrate again that the eigenvalues x\,n of Ji are well separated, hence we can 
expect the accuracy of eigenvectors g] n to stay reasonably close to machine precision. Indeed, the error is below 
120eM for all values of «, as indicated by Figs.[8{d-f), where 6m = 2~ 53 » 1.11 x 10~ 16 denotes the machine epsilon 
in double precision [33, p. 79]. Considering the tridiagonal form of J m with the simple expressions (l98l for the matrix 
elements, its superiority over K,„ in the calculation of eigenvectors is justified. 

4. Concluding remarks 

We have formulated a scalar problem which is equivalent to the concentration problem of tangential vector fields 
within a spherical cap, and enables us to treat it analogously to the concentration problem of scalar functions. Hence a 
construction of a commuting differential operator with a simple spectrum has been made possible. This circumstance, 
at the same time, opens the way for computing concentrated vector fields in a fast and numerically stable way, as 
opposed to the direct method based on the ill-conditioned concentration matrix. 

The reduction of the vector problem to an equivalent scalar one relies on a special combination of vector spherical 
harmonics, which we used as basis functions throughout this paper. With the help of the functions F\ m of Sheppard 
and Torok, for which we derived several novel relations, our mixed vector spherical harmonics can be expressed in 
a simple separable form. Finally, we note that these novel relations of F\ m could facilitate the development of a fast 
vector spherical harmonic transform, too fl34n . 
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Appendix A. Proofs 

Appendix A.l. Proof of orthononnality relation ( Tl9b 

We first substitute definition (TTTb of F\ m into the left-hand side of orthonormality relation ([T9V This yields 



X 



F hn (x)Fi> m (x) Ax = 
1 '(< + 



2 dUi m (x) dUi> m (x) m Ui m (x)Ui> m (x) 
(l-x ) — : : + 



i r , 

+ 1) J_j [ dx dx 

r 1 d 



l-x 2 



dx 



The first term evaluates to [22, p. 754] 
/(/ + _ ' 

while for the second term, we get 
m 



— r 
+ 1) j-i 

we get 
m f 1 d 



1(1+1) J_ 



2 .dUi m (x)dUy m (x) m 2 U Im (x)Up m (x) 
(l-x ) — : : + 



dx 



dx 



l-x 2 



dx = 6 W . 



(x)U Vm (x)] dx 



-[U, m (\)U Vm (\) - Ui, n (-l)Ur m (-l)] = 0, 



i(i+ iy 

where the last equality follows from the expression rtl4t for the values of t// m (±l). Thus we end up with 

l 

Fi m (x)F Vm (x) dx = 5w ■ □ 



Appendix A.2. Proof of recurrence relation ( 1231 ) 

Using expressions (|20T > and (fJTJ of F\ m and recurrence relation (0 of C/; m , we transform the left-hand side (LHS) 
and right-hand side (RHS) separately so that only terms containing Ui,„ and C/;-i im remain. 

First we rewrite the LHS by inserting (l2(Tt : 

-(lx + m)Ui m (x) + (21 + l)£ m t//_i, m (x) 



L1IS = X "7<7TT)J F "" W = 

(lx + m)[l(l + l)x - m] 



1(1 + 1) 



y/l(l+ 1) Vl -x 2 



(2/+l)[/(/+l)x-m] 

!/|.,W + , : -%lmUl-\,m(x). 



[i(i + 1)] 3/2 VT^ 5 [/(/ + 1)] 3/2 VTT? 

After that we proceed to the RHS. We insert (|20T > and (EH . shifted in index / by +1 and — 1, respectively: 

RHS = (i+\, m Fi+\, m (x) + £lmFl-\,m(X) 

-[(I + l)x + m] Ui+i, n (x) + (21 + 3)£ /+ i, m t/; m (x) (lx - m)Ui- Um (x) - (21 - Y)i;i m Ui m (x) 



V(/ +!)(/ + 2) VT 



+ 4n- 



V(/- i)/Vi -x 2 



Next we expand ^/+i >m and using their definition ( I26l l. By straighforward, if lengthy, algebraic calculation, we get 

DUC _ ~l 2 [(l + 1)JC + mfoufiti,«W + (' + 1) 2 (^ - m)ti m Ui- Um (x) + m 2 Ui m (x) 
KHS - — - — . 

[i(i + 1)] 3/2 VT^l? 

We apply recurrence relation (1231 and collect like terms, hence 

-/ 2 x[(7 + l)x + m]+m 2 l 2 [(l + l)x + m] + (I + l) 2 (/x - m) 

RHS = — , -Ulm(x) + ; — , &mUl-l, m (x). 



[/a+i)] 3/2 vr^ 



[id + 1)] 3/2 vr^ 



Taking the difference LHS - RHS, it can be shown by further straightforward algebra that the coefficients of Ui, n and 
Ui-\ m are zero. Hence LHS = RHS. □ 



15 



Appendix A.3. Proof of recurrence relation ( 1241 ) 

In this proof, we follow the same strategy as in the previous proof and rewrite the left-hand side (LHS) first. 
Inserting expression (l20t of F\ m yields 



LHS = (1 - = — — n - 

dx + 1) dx 



-(/x + m)UiJx) + (21 + l)£, m t//_i >m (x) 



vr^ 2 



Performing the differentiation and using (1 - x 2 )^(l - x 2 ) '^ 2 = x(l - x 2 ) ^ 2 , we get 

- (Ix + m)xUi m (x) + (21 + l)£/ m xt/;_i jm (x) - 1(1 - x 2 )Ui m (x) - (Ix + m)(l - x 2 ) 



LHS = 



dUlm(x) 

dx 



+ (21+ 1 )£,„(! -x 2 ). 



dx 



x [/(/ + 1)(1 - x 2 )] 



2x-|-l/2 



Next we insert recurrence relations (fj} and ([?}, shifted in index Z by +1 and - 1 , respectively. After that we collect like 
terms and perform some straightforward algebra to obtain 

(Zx + m)(Z-l)x-Z(l-x 2 )-Z 2 +m 2 rr (2Z+l)(x-m) tj 

LHS = — = -= = C//m(x) + £unUl-\,nA.X). 

VZ(Z + 1) Vl - x 2 VW+ 1) Vl -x 2 

Now we rewrite the right-hand side (RHS). We insert expressions ( T20t and (I2T1) of F; m , the second one shifted in 
index Z by + 1 . 



(in \ 
x--]Fi m (x) + (2l+l)£ lm Fi. 



-l,m(x) 

(w/Z - Zx)[-(Zx + m)Ui m (x) + (21 + \)^ m Ui-i M (x)] (21 + l)[(Zx - m)U,- Un (x) - (21 - l)£ m t/ /m (x)] 

_ ; h * ; ■ 

VZ(Z+ l) Vl -x 2 V<7- l)/ Vl -x 2 

Next we substitute definition (l26T i of £/ m and collect like terms. By straightforward algebra we get 



(lx + m)(l 2 x-m)-(l + l)(Z 2 -m 2 ) (2Z+l)(x-m) 

RHS = Ul m (x) + ; 6mt//-| ,„(x). 



zy^ + 1) Vi -x 2 



VZ(Z+ i) Vl -x 2 ' 



Taking the difference LHS - RHS, the terms containing £/;-i,m cancel. It can be shown that the coefficient of Ui m 
is zero as well, hence LHS = RHS. □ 

Appendix A.4. Proof of Christojfel-Darboux formula ( 1271 ) 

We start from recurrence relation (ETl and multiply both sides by Fi m (x'). Then we take the same recurrence 
relation again, but this time, substitute x' for x and multiply both sides by Fi m (x). In this way, we obtain the following 
two equations: 



m 



1(1+ 1) 
m 



1(1+1) 

Taking their difference and summing over Z yields 



F lm(x)F l m (x') - £l+\, m Fl+\jn(x)Flm(x) + QmF l-\, m (x)F i m (x'), 
F lm(x')F i, n (x) = (l+l, m Fl+ljn( x ')Flm(x) + ClmF l-\j n (x')F t m (x). 



(x - x) Fi m (x)Fi, n (x) = ^j^i + un[Fi +Um (x)F, m (x) - Fi m (x)Fi +Um (x')\ 



1=1 m 



l=C„ 



■ Clm[Fl-Um(x)Fi m (x') - F / ffl (x)F/_i im (x')]}. 



We can see that consecutive terms cancel in the sum on the right-hand side. Moreover, F( m -\ >m = 0, thus only one 
term corresponding to remains: 



(X - x') Fi m (x)F lm (x) = (L+ljn[FL+\,m(x)FLm(x) - F Lln (x)F L+h ,„(x')] . □ 
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'=4, 



Appendix A. 5. Proof of F i m satisfying differential equation (l28T l 
First let us rearrange (EFt and insert F[ m : 



d 



-/(/+ 1) + 



to 2 - 2tox + 1 



1 -x 2 



Flm(x). 



(A.l) 



The left-hand side can be transformed by exploiting recurrence relations (l24l i and (l25l l (the second one shifted in index 
/ by -1) as follows: 



d 

6x 



n 2 .dFi m (x) 
(1 -x )■ 



dx 



d 

dx 



(to \ 
x - -pJFi m (x) + (21 + \)^ m Fi^ m (x) 



-l 2 (l - x 2 )F lm (x) - (l 2 x - m)(l - x 2 )- dF "" ( Vi 



dx 



+ /(2/- lkWl 



dx 



2M-1 



x[/(l-x 2 )] 



= {-/ 2 (1 - x 2 ) J F to (x) - (/ 2 x - m) -l(x - j)F lm (x) + (21 + 1)^. 

l(x - '^Fi- hm (x) ~ W ~ m,nFim(x) } x [1(1 - x 2 )] 



-l,m(x) 



+ 1(21+ 1)6 



2 VI -1 



Collecting like terms yields 



d 

dx 



2 dF ;m (x) 
(1 - x ) 



dx 



-/ 2 (1 - x 2 ) + /(7 2 x - m)(x - m/l 2 ) - 1(21 - l)(2l + \)tf m 



1(1 -x 2 ) 

-(21 + l)(l 2 x - to) + (21 + l)l 2 (x - m/l 2 ) 
1(1 -x 2 ) 



Fi m (x) 



(lmFl-l, m (x). 



As expected, the term containing F;_i >m vanishes. Applying definition d26l l of £/ m and expanding the fraction by I 
yields 



d 

dx 



2 dF[ m (x) 
(1 -x )■ 



dx 



-Z 3 (l - X 2 ) + (/ 2 X - TO) 2 - (I 2 - l)(l 2 - TO 2 ) 



Hi -x 2 ) 



F, m (x). 



By a straightforward, if lengthy, simplification we obtain 



d 

dx 



2 -, dF , m (x) 
(1 -x )■ 



dx 



-1(1 + 1) + 



to 2 — 2tox + 1 



1 -x 2 



Fi m (x). □ 



Appendix A.6. Proof of addition theorem ( 133b 

Upon inserting definition ([T7| > of F; m into the left-hand side of addition theorem Q31> . we obtain 

m--l \m=-l L J m—-l ) 



— t 

+ 1) A 



i(i + 1) 



dt// m (x) 
I 

dx 



Ul, n (x). 



The first term yields (21 + 1)1 2 because of addition theorems ( TToT l. while the second term can be proven to vanish as 
follows: 



dUiJx) 

il, It / i m A U lm(x) = 

1(1 + 1) dx 



— Y 

'+ 1) A 



m=—l 



i(i + 1) 

2 



ZdUiJx) vi dUi m (x) 
to — U tm (x) + 2^ m — Uim(x) 



m=-l 
I 



m=l 



/(/ 



— y 

+ 1) -H 



d£/ f ,_,„(x) dt/ fa ,(x) 
(-to) — Ui- m (x) + m — — — Ui m (x) 



dx 



dx 
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Using symmetry relation (TTZt . we get 



z y m dU^l = <_ y [ ( _ 1)2m( _ w) + m ] 
+ 1) dx /(/+ 1) L J 

m=-/ m=l 



/(/ + 1) ^, dx 

m=-l 



dUlm(x) 

dx 



Ulmix), 



hence 



Therefore, 



2 — V 



dt// m (x) 
dx 



t// m (x) = 0. 



2 

m=—l 



2 21+1 



□ 



Appendix A.7. Proof of integral identity d9"TT > 

Inserting expression ((90} of J7m into both sides of integral identity d9"TT l yields 



it\(x) 



\ ui(x)[J m u 2 (x)]dx = I 

«7cos0 Jcos© 

+ /' 

«7cos 

| [J m ui(x)]u 2 (x)dx = f 

Jcos© Jcos0 UX 

Jcos 



dx 



(COS © - x)(l - X ) 

dx 



L(L + 2)mi(x)«2(x) - (cos® - x) 
(cos©-x)(l -x 2 ) 



2x d«i(x) 



dx 



L(L + 2)u\(x)u2(x) - (cos - x)- 



dx 

m 2 — 2mx + 1 
1 -x 2 

M2(x) dx 

m 2 — 2mx + 1 



Mi(x)i<2(x) 



1 -x 2 



-u\(x)ui(x) 



Next we perform integration by parts on the first term of the right-hand side in both equations: 



f 1 d 

mix)— 

Jcos© ax 



(cos0-x)(l -x 2 ) 



dt<2(x) 



dx 



dx = Mi(x)(cos - x)(l - x 2 ) 



d«2(x) 
dx 



f - 

Jcos0 dx 



«7cos 



dwi(x) 2 dM 2 (x) 

(cos © - x)(l - x ) — : dx, 



dx 



dx 



„ n dMr(x) 

(cos©-x)(l -x 2 )— — 
dx 



T dMl(x) 

K2M dx = (cos - x)(l - x ) ui(x) 



dx 



s© 



- r (cos©- 

Jcos© 



2 dMi(x)dM 2 (x) 

X)(l - X ) ; ; dx. 



dx dx 

The first term on the right-hand side of both equations vanishes and the rest is identical, hence 

- 1 d 

dx 



«7cos 



d 

u\{x) 



dx 



duo(x) 
(cos0-x)(l -x 2 ) " 



dx 



dx 



f 

Jcos© 



r, dMr(x) 

(cos ©-*)(! -x 2 ) 



dx 



Upon inserting (IA.5b into (lA.4at we find that 



«7cos 



u\{x) [J m u 2 (x)]dx 



-f 

Jcos© 



[J' m u l (x)]u 2 (x)dx. □ 



(A.2) 



(A.3) 



dx, (A.4a) 



dx. (A.4b) 



U2(x)dx. (A. 5) 
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Appendix A.8. Proof of identity ( 1931 ) 

First we apply expression ((89} of J7m to the kernel function "7C m (x, x') and use eigenvalue equation d43T > of An , 



J m %„(x,x') = 



(cos - x)A S i m - (1 - x ) L(L + 2)jc 



2^ Fi m (x)F hn (x') 

i=e„, 



- cos J] /(/ + l)F /m (*)F ;m (*') + x + 1) - L(L + 2)]F to (x)F k (/) 



i=4, 



-d-- 2 )£ 



dFi m (x) 
dx 



Flm(x'). 



Likewise, we also apply $' m to 'Kix, x') and subtract the resulting equation from the previous one, yielding 

L 

or* - j^u, *o = u - *o 2 + 1} - L(L + ^i^f*^ 1 ') 



dF,,„(x) 



L 



Ax 



Flm(x') 



+ Y,F, m (x)(l-x' 2 ) 



l2 dFi m (x') 



l=tm 



Ax' 



Using recurrence relation (1251 on the terms containing the derivatives of F/ m and performing some straightforward 
algebra, we get 

L 

(J m ~ J'JKJx, x') = (jc - x') Yjl 1 ~(L+ l) 2 ]F lm (x)F lm (x') 

i=(„ 

L 

+ 2^(2/ + lXi+i, m [F M , m (x)Fi m (x') - Fi m (x)Fi +hm (x')]. 



Applying Christoffel-Darboux formula (l27t to the second term on the right-hand side yields 

L 

(J m ~ J' m YK m (x, x') =(x- x') J] F, m (x)F lm (x')[l 2 - (L + l) 2 ] 

L / 

+ (x - x') Yp-l + 1) Yj F i'm{x)F Vm {x'). 
i=t m i'=e„, 

In the last term of the right-hand side, the summation can be interchanged as 

l i L L 

Yp-l +1)2 F,, m (x)F,, m (x') = ^ F Vm {x)F Vm {x') Yp-l + 1). 



(A.6) 



/=/' 



Relabeling the sums on the right-hand side of this expression, so that / becomes ¥ and vice versa, and inserting the 
resulting expression into the right-hand side of (1A.6K we obtain 

L r L 

(J m - J' m W m {x, x) = {x- x') J] F lm (x)F lm (x) l 2 -(L+ l) 2 + J](2/' + 1) 



1=1, „ 



l'=l 



Since YTi>=t + 1) = (L + 1) - I , the right-hand side vanishes. Hence 

J m K m {x,x')=J' m K m {x,x'). □ 
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Appendix A.9. Proof of the matrix elements d98t of Sf m 

We start by inserting expression (f89t of J~ m into the integral expression (|96l for the matrix elements and use the 
eigenvalue equation d45b of An,,,: 



(cos - x)Aa„, - (1 - x ) L(L + 2)x 

ax 



F Vm (x) Ax 
1 



= -/'(/' + 1) cos J Fi m (x)F,, m (x) dx + [/'(/' + 1) - L(L + 2)] J" xFi m (x)F Vm (x) Ax 



;_i dx 



(A.7) 



The first integral is equal to 6w because of orthonormality relation ([T9V The remaining two can be evaluated by using 
recurrence relations (IZTt and (124} and orthonormality relation ( fT9T >: 



x: 



x: 



xFi m (x)F Vm {x)&x = £r m 8m-i + £> + i ; „A/+i + yy, ~ yj fa' 



F lm (x)(l-x z ) 



2 &F Vm (x) 



dx 



m 

dx = (/' + \)£i> m 8i t i>-\ - V ' £i>+i, m 6i,i>+i + — - 6 iv 



Thus for (IA.71 >. we get 



J m ,iv = \ -I' (I' + 1) cos + m 



1 - 



L(L + 2) + 1 



I'd' + 1) 

+ ^ + i,m[/'(/'+2)-L(L + 2)]<J r/+1 . 



5//' + &m [(/' - + 1) - UL + 2)] 



Because of the Kronecker deltas, this expression is non-zero for index pairs (I, I), (I + 1, /) and (/, I + 1) only. The 
corresponding matrix elements are 



m,ll 



-1(1 + l)cos0 + m 



1 



L(L + 2) + 1 



/(/ + 1) 

Jm,l,l+l = Jm,l+l,l = f/ti#[l(l + 2) - L(L + 2)], 



hence J m is real, symmetric and tridiagonal. □ 
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Figure 2: Sketch of the unit sphere CI showing Cartesian unit vectors x, y, z, and the polar and azimuthal unit vectors and 0, respectively. 
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Figure 5: Rank-ordered eigenvalue spectrum including the eigenvalues rj m of all K,„, — L < m < L for (a) = 30°, (b) © = 60°, (c) © = 90°, and 
L = 18. The vertical gridlines mark the corresponding Shannon numbers N of {85). 
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Figure 6: Four scalar eigenfunctions G m „(cosff), n = 1,3,5,7 of each order -2 < m < 2. The maximal degree is L = 18 and 6) = 60°. The 
black and gray curves mark contributions of G,„„ to the interior of the spherical cap (0 < 8 < 60°) and the rest of the sphere (60° < 9 < 180°), 
respectively. Labels show the eigenvalues rj nm which express the quality of concentration within C. 
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Figure 7: Rank-ordered eigenvalue spectrum including the eigenvalues Xm of all J,„, -L < m < L for (a) = 30°, (b) = 60°, (c) = 90°, and 
L = 18. The vertical gridlines mark the corresponding Shannon numbers N = 24, 90, 180 (see (85}). 
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Figure 8: (a-c) Eigenvalue gap i ll 00) for Ki (circle) and Ji (triangle) for = 30°, 60°, 90°, respectively, (d-f) Error d 1 1 1 of the eigenvectors of 
Ki (circle) and Ji (triangle) for = 30°, 60°, 90°, respectively. The number £m denotes the machine epsilon in double precision. The maximal 
degree is L = 18 and the vertical gridlines mark the partial Shannon numbers Ni of (84). 
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